Event-by-Event Fluctuations in Particle Multiplicities and Transverse Energy Produced in 

158 -A GeV Pb+Pb collisions 
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Event-by-event fluctuations in the multiplicities of charged particles and photons, and the total trans- 
verse energy in 158 -A GeV Pb+Pb collisions are studied for a wide range of centralities. For narrow 
centrality bins the multiplicity and transverse energy distributions are found to be near perfect Gaus- 
sians. The effect of detector acceptance on the multiplicity fluctuations has been studied and demon- 
strated to follow statistical considerations. The centrality dependence of the charged particle multi- 
plicity fluctuations in the measured data has been found to agree reasonably well with those obtained 
from a participant model. However, for photons the multiplicity fluctuations have been found to be 
lower compared to those obtained from a participant model. The multiplicity and transverse energy 
fluctuations have also been compared to those obtained from the VENUS event generator. 
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I. INTRODUCTION 

Fluctuations in physical observables in heavy ion col- 
lisions have been a topic of interest for some years as 
they may provide important signals regarding the forma- 
tion of Quark-Gluon Plasma (QGP) and help to address 
the question of thermalization [|J. With the large num- 
ber of particles produced in heavy ion collisions at SPS 
and RHIC energies ^,||], it has now become feasible to 
study fluctuations on an event-by-event basis. Recently, 
several new methods have been proposed for the study of 
event-by-event fluctuations in various global observables 
to probe the nature of the QCD phase transition [Q-^]. 
In a thermodynamical picture of a strongly interacting 
system formed in the collision, the fluctuations in parti- 
cle multiplicities, mean transverse momenta ((pt)), and 
other global observables, are related to the fundamental 
properties of the system, such as the specific heat, chem- 
ical potential, and matter compressibility. These, in turn, 
lead towards understanding the critical fluctuations at the 
QCD phase boundary. The existence of a tri-critical point 
at the QCD phase transition which has lately been a 
topic of intense discussion, has been predicted to be asso- 
ciated with large event-by-event fluctuations in the above 
observables. 

In a first order phase transition scenario, it is be- 
lieved that supercooling might lead to density fluctua- 
tions resulting in droplet formation and hot spots [||]. 
These might lead to rapidity fluctuations in the form of 
spikes and gaps in the rapidity distribution. The study of 
event-by-event fluctuations in the number of photons to 
charged particles has also been proposed as a means to 
search for production of Disoriented Chiral Condensates 
(DCC) HQ. 

In nucleus - nucleus collisions the transverse energy, 
Et, is an extensive global variable jrT|-|l3|] which pro- 
vides a direct measure of the violence of an interaction. 
Et is produced by redirection of the longitudinal energy 
into transverse motion through interactions in which the 
interacting particles undergo multiple scatterings and ap- 
proach thermalization. Et is also an indicator of the en- 
ergy density achieved in the collision. Since the energy 
density is directly related to the QGP phase transition, it 
is extremely important to study Et and fluctuations in 
Et- Moreover it is interesting to compare the fluctua- 
tions of Et to those observed in the particle multiplici- 
ties. 

Much theoretical interest has been directed toward the 
subject of event-by-event fluctuations, motivated by the 
near perfect Gaussian distributions of (pt) and parti- 
cle ratios Jh] ] measured at the SPS. For these Gaussian 
distributions, the variance or the width of the distribu- 



tions contain information about the reaction mechanism 
as well as the nuclear geometry [H,|l5- 17 1. 



The relative fluctuation (u>x) m an observable X can 
be expressed as: 



{xy 



(i) 



where a x is the variance of the distribution and (X) de- 
notes the mean value. 

The value of ujx which can be extracted from experi- 
mental data has contributions which originate both from 
trivial statistical effects as well as dynamical sources. To 
extract the dynamical part associated with new physics 
from the observed fluctuations, one has to understand the 
contributions from statistical and other known sources. 
Examples of known sources of fluctuations contributing 
to the observed experimental value of u>x include finite 
particle multiplicity, effect of limited acceptance of the 
detectors, impact parameter fluctuations, fluctuations in 
the number of primary collisions, effects of re-scattering 
of secondaries, resonance decays, and Bose-Einstein cor- 
relations. These sources of fluctuations, along with esti- 
mates of the lox contributions for each have been dis- 
cussed by Stephanov et al. [Qj and by Heiselberg et al 

In nucleus-nucleus (AA) collisions relative fluctua- 
tions in global observables have been found to be smaller 
compared to those in pp collisions. It is suggested that 
thermal equilibration in AA collisions makes the fluc- 
tuations small. However, the origin of fluctuations and 
hence the physical information content are quite differ- 
ent in pp and AA collisions. While in pp collisions 
one hopes to extract quantum mechanical information 
about the initial state from the event-by-event fluctua- 
tions in the final state, in heavy-ion collisions equilibra- 
tion makes it difficult to achieve this goal, instead the 
basic aim here has been to relate the event-by-event fluc- 
tuations of the final state with the thermodynamic prop- 
erties at freeze-out. 

In this article we present fluctuations in the multiplic- 
ities of both charged particles and photons, and in the to- 
tal transverse energy, over a large range of centralities as 
measured in the WA98 experiment at the CERN SPS. A 
major interest has been to search for fluctuations which 
have a new physics origin, such as those associated with 
QCD phase transition or from the formation of a DCC. 

We compare the fluctuations observed in the experi- 
mental data for varying centrality conditions and rapidity 
intervals to those obtained from different models. In the 
next section the WA98 experimental set up is described. 
In section 3 the criteria for the centrality selection appro- 
priate for fluctuation studies are discussed. Multiplicity 
fluctuations of photons and charged particles and the ef- 
fect of acceptance are presented in section 4. In section 
5, we estimate the multiplicity fluctuations in a partic- 
ipant model and compare to those obtained from data. 
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Section 6 deals with transverse energy fluctuations. A 
final discussion and summary is presented in section 7. 

II. EXPERIMENT AND DATA ANALYSIS 

In the WA98 experiment at CERN 0, the main em- 
phasis has been on high precision and simultaneous de- 
tection of both hadrons and photons. The experimental 
setup consisted of large acceptance hadron and photon 
spectrometers, detectors for charged particle and photon 
multiplicity measurements, and calorimeters for trans- 
verse and forward energy measurements. The exper- 
iment took data with 158-^4 GeV Pb beams from the 
CERN SPS in 1994, 1995, and 1996. The results pre- 
sented here are from the Pb run in 1996 taken with the 
magnet (Goliath) turned off. The analysis makes use 
of the data taken with the photon multiplicity detector 
(PMD), the silicon pad multiplicity detector (SPMD), the 
mid-rapidity calorimeter (MIRAC), and the zero degree 
calorimeter (ZDC). 

The circular Silicon Pad Multiplicity Detector 
(SPMD), used for measuring charged particle multiplic- 
ity, was located 32.8 cm from the target. It had full az- 
imuthal coverage in the region 2.35 < i] < 3.75. The 
detector had four overlapping quadrants, each fabricated 
from a single 300 fim thick silicon wafer. The active 
area of each quadrant was divided into 1012 pads form- 
ing 46 azimuthal wedges and 22 radial bins with pad 
size increasing with radius to provide a uniform pseudo- 
rapidity coverage. The intrinsic efficiency of the detector 
was better than 99 %. During the data taking, 95 % of the 
pads worked properly. It was nearly transparent to high 
energy photons, since only about 0.2 % are expected to 
interact in the silicon. Details of the characteristics of the 
SPMD can be found in Ref. [^9@]. 

The photon multiplicity was measured using the 
preshower photon multiplicity detector (PMD) placed at 
a distance of 21.5 meters from the target. The detector 
consisted of 3 radiation length (Xq) thick lead converter 
plates placed in front of an array of square scintillator 
pads of four different sizes, varying from 15 mmx 15 mm 
to 25 mmx 25 mm, placed in 28 box modules. Each box 
module had a matrix of 38 x 50 pads which were read out 
using one image intensifier + CCD camera system. De- 
tails of the design and characteristics of the PMD may be 
found in Ref. [ pTp^ ]. The results presented here make 
use of the data from the central 22 box modules covering 
the pseudo-rapidity range 2.9 < rj < 4.2. The clusters 
of hit pads, having total ADC content above a hadron re- 
jection threshold were identified as photon-like. Detailed 
simulations showed that the photon counting efficiencies 
for the central to peripheral cases varied from 68% to 
73%. The purity of the photon sample, A r 7 _i;] <0 , in the 
two cases varied from 65% to 54%. 

The transverse energy was measured with the MIRAC 
calorimeter [|23|] placed at 24.7 meters downstream from 



the target. It consisted of 30 stacks, each divided ver- 
tically into 6 towers, each of size 20 x 20 cm 2 , and 
segmented longitudinally into an electro-magnetic (EM) 
and a hadronic section. The depth of an EM section 
was 15.6Xo (equivalent to 51% of an interaction length) 
which ensured almost complete containment of the elec- 
tromagnetic energy (97.4% and 91.0% containment cal- 
culated for 1 GeV and 30 GeV photons, respectively). 
The MIRAC was used to measure both the transverse 
electromagnetic (E^ 71 ) and hadronic (E^ ad ) energies 
in the interval 3.5 < 77 < 5.5 with a resolution of 
17.9%/VE and 46.1%/ y/E, (E in GeV), respectively. 
The Et provides a measure of the centrality of the re- 
action. Events with large Et correspond to very central 
reactions with small impact parameter and vice versa. 

The Zero Degree Calorimeter (ZDC) measured the to- 
tal forward energy, Ep, at 6 < 0.3° with a resolution 
of 80%/a/E + 1-5%, with E expressed in GeV. E F pro- 
vides complementary information on the centrality with 
low E-p energy deposit corresponding to small impact pa- 
rameter collisions. 

For the results to be presented below, the following 
sources of errors have been included in the systematic 
error estimates: 

Errors in N 1 -\\\ ic : 

(a) The major source of error in N^^\^ c is due to the 
effect of clustering of the pad signals. This error is de- 
termined from the simulation by comparing the number 
of known tracks on the PMD with the total number of 
photon-like clusters. The result is that the number of 
clusters exceeds the number of tracks by 3% in the case 
of peripheral events and by 7% for high multiplicity cen- 
tral events [53]. 

(b) The uncertainty in the ADC value of the hadron rejec- 
tion threshold in the PMD leads to an error in the estima- 
tion of A^-iikc clusters. The hadron rejection threshold 
has been set at three times the MIP (minimum ionizing 
particle) peak. The value of MIP peak was changed by 
10% of the peak value (3 ADC) in order to estimate the 
systematic error. The error in iV 7 _iikc value is 2.5% fl22|]. 

(c) The error because of the variation in scintillator pad- 
to-pad gains is found to be less than 1%. 

The combined systematic error on -/V 7 _iik c is asym- 
metric and centrality dependent. The errors are -3.2% 
and +3.4% for peripheral collisions and -7.1% and 
+3.0% for central collisions. The errors on N y , obtained 
after correcting for photon counting efficiency and purity 
of photon-like sample, will be discussed in section V. 

Errors in N c h : 
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The uncertainty in the obtained from SPMD has 
been discussed in detail in Ref. Jl^]. The total error has 
been estimated to be about 4%. 

Errors in centrality selection through Et and Ep : 

The centrality of the interaction is determined by the 
total transverse energy (Et) measured in the MIRAC. 
The finite resolution in the measurement of Et con- 
tributes to the systematic error. For the analysis of fluc- 
tuation in Et which uses MIRAC data directly, the cen- 
trality is determined by the forward energy, E-p. The fi- 
nite resolution in the measurement E-p contributes to the 
systematic error in Et fluctuations. These errors are cen- 
trality dependent. 

Fitting errors : 

The fitting errors associated with the determination of 
the fit parameters of the multiplicity and transverse en- 
ergy distributions also contribute to the final systematic 
error in both the photons and charged particles and trans- 
verse energy respectively. The maximum contribution of 
this error to the fluctuation was found to be 2%. 



III. CENTRALITY SELECTION FOR FLUCTUATION 
STUDIES 

The centrality of the interaction was determined by the 
total transverse energy measured in the MIRAC. For the 
part of the analysis where transverse energy data are used 
for the fluctuation studies, the centrality was determined 
instead by the forward energy, E-p, as measured in the 
ZDC. The centralities are expressed as fractions of the 
minimum bias cross section as a function of the mea- 
sured total transverse energy using MIRAC, or total for- 
ward energy using the ZDC. Fig. |l](a) and (b) show the 
minimum bias distributions of Et and Ep, respectively. 
The arrows in the figures indicate the values of Et and 
E F for the top (most central) 1%, 2%, 5%, and 10% 
of the minimum bias cross section. Predictions from 



VENUS 4.12 [24 1 are also shown as solid histograms. 



This will be discussed in a later section. 

The anti-correlation of Et and Ep is shown in[j](c). It 
illustrates that either Et or Ep can be used nearly equiv- 
alently to define the centrality of the reaction. 

Fig. |(a) and gb) show the minimum bias distribu- 
tions for 7 — like clusters and charged particles, respec- 
tively, for the full acceptances of the two detectors (PMD 
and SPMD). The multiplicity distributions correspond- 
ing to the centrality cuts using the total Et for the top 
1%, 2%, and 5% of the minimum bias cross section are 
also shown in the same figures. These distributions have 
been fitted to Gaussians. The extracted fit parameters are 
used for the analysis of the fluctuations. 



Fig. [}] shows the variation of the mean (n), stan- 
dard deviation (a), and chi-square per degree of freedom 
(X 2 /ndf) of the photon and charged particle multiplicity 
distributions for different centrality bins. Here the cen- 
trality class is chosen with increasing width, as — 1%, 
0-2%, 0-3%,- • -,0-10%. As expected, the mean value 
decreases and the sigma increases as we make broader 
centrality selection to include more of the cross section. 
From the x 2 /ndf values, one observes that the distribu- 
tions increasingly deviate from Gaussians with increas- 
ing width in the centrality selection. For a centrality se- 
lection width of greater than 5%, the x 2 /ndf rises above 
2. The variation of /i and a indicate that the extracted rel- 
ative fluctuation (wj) will grow with the increase in the 
width of the centrality selection interval. This indicates 
that the impact parameter fluctuations will dominate as 
the centrality selection is broadened. From this we con- 
clude that the centrality selections should be made with 
as narrow bins in Et as possible, such that the multi- 
plicity distributions are good Gaussians and the impact 
parameter fluctuations are minimized. With this in mind 
we have used centrality selection bins of 2% widths in 
cross section, taken as, - 2%, 2 - 4%, • ■ •, 62 - 64%. 

Fig. |] shows the variation of u, a, and x 2 /n<i/ of 
the photon and charged particle multiplicity distributions 
within the full acceptance of the detectors with these nar- 
row bins in centrality. The data presented here cover the 
region from central (top 2% of the minimum bias cross 
sections) to peripheral collisions (up to 65% of the min- 
imum bias cross section where the average number of 
participants is 26). It is seen that both the a and a val- 
ues decrease towards peripheral collisions. The x 2 /ndf 
values are mostly in the region between 1.0 and 2.0 over 
the entire range of centralities considered. This suggests 
that narrow cross section slices in the Et or Ep distribu- 
tions are necessary to study the multiplicity fluctuations 
and minimize the influence from impact parameter fluc- 
tuations. 



IV. MULTIPLICITY FLUCTUATIONS AND THE 
EFFECT OF ACCEPTANCE 

The relative fluctuations in multiplicity for 7 — like 
clusters and charged particles have been calculated using 
the n and a values from Fig. |] and Eq. 1. These val- 
ues are shown in Fig. || as functions of centrality, for 
full coverage of PMD (2.9 < n < 4.2) and SPMD 
(2.35 < n < 3.75), respectively. The errors shown in 
the figures are systematic errors, the sources of which 
has been already discussed in previous section. For both 
7 — like clusters and charged particles the relative fluc- 
tuations are seen to increase in going from central to pe- 
ripheral collisions. However, for charged particles the 
increase is much stronger. 

In order to make a direct comparison of the fluctua- 
tions of photons and charged particles, the multiplicities 
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should be studied in the region of common coverage of 
the detectors in terms of both pseudo-rapidity and az- 
imuthal angle. In a later section we will compare the 
results obtained from data with those from model calcu- 
lations for the common coverage. The region of common 
coverage of the two detectors in the WA98 experiment 
was 0.85 units in 77 (2.9 < 77 < 3.75). The general trend 
of the variation of the Gaussian fit parameters, with the 
reduced number of particles, for the common coverage is 
found to be similar to that obtained with full coverage for 
each detector. Fit parameters /i, er, and /ndf, as ob- 
tained for centrality bins of 2% in width, for the common 
coverage of the two detectors are shown in Fig. ^. As was 
the case for the larger acceptance, the ji and a values are 
seen to decrease towards more peripheral event selection. 
The x 2 /ndf values are reasonable. Using the above val- 
ues of the Gaussian parameters together with Eq. 1, the 
relative multiplicity fluctuations were calculated and are 
shown in Fig. ^. The error shown include the fit errors 
as well as the other systematic errors discussed earlier. 
The relative fluctuations for 7 — like clusters is seen to 
be rather constant over the full centrality range with an 
average value of 2.2 ± 0.21. In comparison, the relative 
fluctuations for charged particles is 1.56 ± 0.13 at the 
most central bin (0 — 2%) increasing to 2.8 ± 0.16 for the 
least central bin (62-64%). 

Following this discussion of the fluctuations in the 
multiplicity of photons and charged particles for the full 
acceptance regions and for the regions of common cov- 
erage of the photon and charged particle detectors, we 
now analyze the effect of detector acceptance on the ob- 
served fluctuations in more detail. For this we have taken 
two different 77 coverage regions for each detector. For 
the PMD the 77 ranges chosen are 3.0 < fj < 4.0 and 
3.25 < ij < 3.75 (with full <f> coverage). The resulting 
lo values for the two cases are shown in Fig. ^|(a). Qual- 
itatively, the variation of the fluctuations with centrality 
is similar for both coverages, but the magnitude of the 
relative fluctuations is lower for smaller 77 coverage. 

For the SPMD, the fluctuations were calculated for the 
rapidity intervals 2.35 < 77 < 3.35 and 2.65 < fj < 
3.15. These bins have width of one and one half unit in 77 
around mid-rapidity. The results are shown in Fig. [|(b). 
Qualitatively, the results are again similar to each other 
with the magnitude of the relative fluctuations decreasing 
as the coverage in 77 is decreased. 

The decrease in the relative fluctuations as the accep- 
tance is decreased can be understood in terms of a simple 
statistical picture [p5[]. Assume that there are m parti- 
cles produced in the collision out of which n particles are 
accepted randomly into the detector acceptance. In this 
case, the distribution of n, will follow a binomial distri- 
bution with mean mf and variance mf(l — /) where / 
is the fraction of particles accepted. Therefore the fluctu- 
ations in the number of particles accepted for a fixed (ra) 
number of particles produced is (1 — /). In principle m 
can have an arbitrary distribution as given by P(m) with 



known first and second moments. The fluctuations in the 
n accepted particles out of the m particles produced is 
then given as: 

w„ = 1 - / + /w m (2) 

Thus considering the fluctuations in one unit of rj as 
ui m we can calculate the expected fluctuations for one 
half unit of 77 using the above equation. Here / corre- 
sponds to the ratio of the total number of particles ac- 
cepted in one half unit of 77 coverage to that accepted 
over one unit in 77. For the acceptance regions used, the 
average value of / for photons is about 0.52 and that of 
charged particles is about 0.54. Using these values in 
Eq. 2 we can calculate the expected fluctuations for half 
unit of rj coverage from the results for one unit of 77 cov- 
erage, under the assumption of a binomial sampling. As 
shown in Fig. ^| the empirical calculations almost exactly 
reproduce the observed result in the narrower acceptance 
window for charged particles and agrees reasonably well 
within the quoted errors for photons. 

V. ESTIMATION OF FLUCTUATIONS IN A 
PARTICIPANT MODEL 

In a picture where the nucleus-nucleus collision is 
thought of as the sum of contributions from many sources 
created in the early stage of the interaction, the variance 
of the distribution of any observable has contributions 
from: 

• the fluctuations in the number of sources, largely 
due to different impact parameters. Even if the im- 
pact parameter window is narrowed, density fluc- 
tuations within the nucleus will make this contri- 
bution non-zero. 

• the fluctuations in the number of particles pro- 
duced by each source. Quantum fluctuations in the 
NN cross section can lead to such effects. 

• the fluctuations due to any dynamical process or 
critical behavior in the evolution of the system. 

The contribution from the first two effects leads to 
fluctuations in the number of participant nucleons which 
may be related to the initial size of the interacting system 
before it thermalizes. Resonance decays have also been 
shown to increase the multiplicity fluctuations by a large 
factor [@|5||. 

Following a simple participant model [P, p3| , p^|j2^ , p7] ] , 
the particle multiplicity (of photons or charged particles), 
N, may be expressed as : 

AW 

i=l 
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where N par t is the number of participants and n, is the 
number of particles produced in the detector acceptance 
by the i th participant. On an average, the mean value 
of m is the ratio of the average multiplicity in the detec- 
tor coverage to the average number of participants, i.e., 
(n) = (JV)/(JV part ). 

Thus the fluctuations in N will have contributions due 
to fluctuations in N paYt N part ) and also due to the fluc- 
tuations in the number particles produced per participant 
(u> n ). Again, the fluctuations of n given as u> n , will have 
a strong dependence on the acceptance of the detector. In 
the absence of correlations between the n/s, the multi- 
plicity fluctuations, ujn, can be expressed as 

lo n = oj n + (n)wjv part (4) 

Comparison of data with the results of such model cal- 
culations might reveal the extent to which the principle 
of superposition of nucleon-nucleon (NN) interactions is 
valid in the case of heavy ion collisions. The participant 
model is expected to hold reasonably well for peripheral 
collisions where there are only few NN collisions, while 
for central collisions the particle production gets affected 
by NN scattering, rescatterings between produced parti- 
cles, energy degradation, and other effects. Next we dis- 
cuss the calculation of each of the terms in Eq. 4. 

A. Calculation 0/Wjv part 

The impact parameter fluctuations are reflected in the 
fluctuations in the number of participants. We have es- 
timated this contribution using the VENUS 4.12 event 
generator with default setting. A set of 100K minimum 
bias Pb+Pb events at 158-A GeV was generated for cal- 
culation of the number of participants. To match the cen- 
trality selection of the reaction in simulation to that in 
data, we have carried out a fast simulation in which Et 
from VENUS was calculated within MIRAC coverage 
taking the resolution factors for the hadronic and electro- 
magnetic energies of MIRAC into account. The corre- 
sponding Et distributions for VENUS are shown as the 
solid curve in Fig. [j](a). It is seen that the agreement with 
data is quite reasonable. 

The distributions of N part for the same narrow (2%) 
bins of centrality, as discussed above for the data, are 
well described by Gaussian distributions. Fig. ^| shows 
the variation of /i, a, x 2 /ndf, and relative fluctuation 
wjv p „ t calculated from the fit parameters with the 2% 
bins in centrality. One can see that the relative fluctuation 
in the number of participants, u>N part < is around 1. The 
statistical errors are small and are within the size of the 
symbols. 

The systematic errors shown in the figures, have con- 
tribution from the following sources, which have been 
added in quadrature. 



• Nucleon density distribution: In order to estimate 
the error due to this we have calculated the number 
of participants from VENUS (as shown in the fig- 
ure) and those from FRITIOF. The difference for 
each centrality bin was considered as representa- 
tive of the error [j28|], 

• Finite resolution of Et' Systematic errors due to 
this were calculated by varying the centrality as 
per the MIRAC resolution [g2|]. 

• Fitting errors: Errors associated with the determi- 
nation of the fit parameters of the Gaussian distri- 
butions also contribute to the final systematic error 
in the number of participants. 

The quantity (n) is equal to the ratio of the mean 
charged particle (or photon) multiplicity for a given ac- 
ceptance to the mean number of participants for the 
same centrality bin. Thus the contribution from the term 
(ft)wjv part to the total fluctuations (Eq. 4) can be easily 
obtained. 



B. Calculation ofu„ 

This term gives the fluctuations in the number of par- 
ticles produced per participant. It has a strong depen- 
dence on acceptance as given earlier in Eq. 2 and shown 
in Fig. ||. To calculate uj n as per Eq. 2 we next obtain 
the terms / and u> m . The quantity / is the ratio of the 
number of particles per participant accepted within the 
acceptance of the detector ((n)) to the total number of 
particles produced per participant ((m)). The value of 
(n) for each centrality bin and for a given acceptance can 
be calculated as discussed in the preceding section. To 
obtain the value of (to) we make use of the data existing 
in the literature for nucleon-nucleon (NN) collisions. As 
discussed in Refs. [ p9p0[ ] the mean number of charged 
particles and photons produced in nucleon-nucleon col- 
lisions can be parametrized as a function of cms energies 
(\fs from 2 GeV to 500 GeV ) in the following manner : 

(N ch ) NN = -4.7(±1.0) + 5.2(±0.8) S °- 145 < ±0 - 01 ) (5) 

(N 7 ) NN = -9.9(±2.1) + 8.5(±1.9) S °- 113 ( ±a015 > (6) 

For the 158 A GeV SPS energy discussed here this pa- 
rameterization gives the average charged particle multi- 
plicity to be 7.2 with the corresponding number for pho- 
tons being 6.3. Thus the average charged particle and 
photon multiplicities per participant are 3.6 and 3.15 re- 
spectively. 

In addition, a for the charged particle multiplicity 
in nucleon-nucleon collisions shows a linear depen- 
dence with the average charged particle multiplicity as 
0.576((^ JV ) - 1), as given in Ref. @. 
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This can be used to calculate u) m , which is given as: 

For charged particles at SPS energies this gives a value 
of oj m = 1.8. However, for photons this number is not 
known since there is no similar parameterization. In the 
absence of such a parameterization of a for photons we 
will also assume that u> m = 1.8 for the photon multiplic- 
ity. Fluctuations of photons, in principle, are expected to 
be similar to those for charged particles. This is because 
the majority of photons come from decay of ir° and while 
the majority of charged particles are charged pions (tt^). 

From the values of (n), (to), and u> m for a given ac- 
ceptance and centrality, the term u> n can then be calcu- 
lated. 



C. Comparison of data to model calculations 

We first compare the experimental results of multi- 
plicity fluctuations to those of the calculations using the 
participant model in the common coverage of PMD and 
SPMD. Fig. p~0| shows a comparison of the fluctuation in 
charged particles from data to that obtained from the cal- 
culations using the model described above. The results 
are plotted as a function of the number of participants. 
The horizontal errors on the number of participants are 
shown only on the data points. The error on ui calculated 
in the model is mainly due to the error on the mean num- 
ber of charged particles in nucleon-nucleon interactions, 
the error in the number of participants calculated, and the 
uncertainty in the simulation of the calculated transverse 
energy. For clarity of presentation we have given results 
corresponding to alternate 2% centrality bins, i.e — 2%, 
2 - 4%, • • •, 62 - 64%. The results from VENUS are also 
shown in form of solid line in Fig. and are found to 
remain almost constant over the entire centrality range. 
Charged particle fluctuations determined from data and 
the participant model decrease in going from peripheral 
to central collisions, although the dependence on central- 
ity is weaker for the model calculation. 

Fig, pi] shows fluctuations in the 7-like clusters as well 
as N~f after the correction. The results, plotted as a func- 
tion of the number of participants, are compared to those 
of the participant model calculations for photons and re- 
sults from VENUS. Using the estimated values of effi- 
ciency (e 7 ) and purity the number of photons in an 
event is calculated by using the relation: 

^V 7 = — ^Vnkc (8) 
e 7 

The photon counting efficiency in PMD varies from 68% 
to 73% for central to peripheral collisions. The purity of 
the measured photon sample varies from 65% to 54% for 
central to peripheral collisions. 



The systematic errors associated with iV 7 _iik c have al- 
ready been discussed in section II. The additional errors 
in the conversion from iV^_uite to N 7 are mainly due 
to errors in estimation of photon counting efficiency and 
purity. These sources of these errors are given below: 

• Event-by-event variation in photon counting ef- 
ficiency (e 7 ) and purity of photon sample (f p ). 
These have been found to vary from 3% to 6% for 
central to peripheral collisions. 

• The purity factor depends on the ratio of the 
number of photons and charged particles within 
the PMD coverage. The systematic error asso- 
ciated with this ratio has been studied by us- 
ing the FRITIOF event generator in addition 
to VENUS. The average photon multiplicity es- 
timated by using FRITIOF is found to be higher 
by about 4% in peripheral and by 1% in central 
collisions, compared to the values obtained using 
VENUS. 

• The photon counting efficiency determined in the 
present case relies on the energy spectra of pho- 
tons as given by the VENUS event generator. As 
the conversion probability for low energy photons 
falls sharply p2| ] with decreasing energy below 
500 MeV, the estimate of e 7 may be affected if 
the energy spectra in the actual case is different. 
Photon energy spectra have been measured by the 
WA98 lead glass calorimeter. By extrapolating 
these measured spectra to the PMD acceptance we 
have estimated the photon counting efficiencies for 
different eta bins and centralities. These results 
turn out to be lower compared to those obtained 
from VENUS by 2-9% for central events and 3- 
13% for peripheral events, the smaller value cor- 
responds to larger pseudo-rapidity region of the 
PMD acceptance. The average difference in ef- 
ficiencies within the PMD acceptance are 6% for 
central and 9% for peripheral collisions. These dif- 
ferences add to the systematic errors on the photon 
counting efficiency. 

The total systematic error on the multiplicity of photons 
(iV 7 ) are -6.7% and +12.5% for peripheral collisions and 
-8.0% and +9.0% for central collisions. 

The fluctuations in the number of photons have been 
estimated from the fluctuations in the number of photon- 
like clusters by using Eq. 8: 

lu 1 — — oj 7 _iiko (y) 
e 7 

These results are shown in Fig. O. It is observed that 
the relative fluctuations of photons from the data are in 
reasonable agreement with those obtained from VENUS. 
However, the results for photons from the participant 
model are somewhat higher than those from the exper- 
imental data. 
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VI. TRANSVERSE ENERGY FLUCTUATIONS 

Relativistic nuclear collisions are often described 
within the participant-spectator picture in which nuclei 
are spheres that collide with a definite impact parameter. 
The overlapping volume which participate in the reac- 
tion are violently disrupted while the remaining spectator 
volumes shear off and suffer comparatively mild excita- 
tions. The magnitude of the Et produced depends on the 
bombarding energy and the participant volume or equiv- 
alently the number of participating nucleons. The cross 
section for a specific value of Et production depends to 
a large extent on the geometric probability of a given im- 
pact parameter. Therefore impact parameter fluctuations 
are expected to lead to fluctuations in Et ■ 

Corroboration of the participant-spectator picture 
comes from the strong anti-correlation of Et with the 
energy observed in the zero degree calorimeter, Ep as 
shown in the Fig. |j](b). The smaller the impact parame- 
ter, the larger is the participant volume and Et, but the 
smaller is the spectator volume and E-p. Et also cor- 
relates strongly with the produced particle multiplicity. 
The da/dNch and da/dN^^u^ spectra have virtually 
the same shape as da/dET (Fig. ||(a) and Fig. ||(b)). 

The study of the average total Et as measured in the 
WA98 experiment and its scaling behavior with the num- 
ber of participants have been discussed earlier in detail in 
Ref. p^]. Here we concentrate on the second moment, 
and study the fluctuations in Et as was done for N^-n^c 
and N c h- For this analysis we have again taken 2% width 
bins in centrality using the forward energy as measured 
by the ZDC. Due to the poorer resolution in centrality se- 
lection of the ZDC for peripheral collisions, we present 
the results only up to the 50% centrality class. 

The Et distribution for the top 2% of the minimum 
bias cross section is shown in Fig. |l2| The solid curve 
shows a Gasussian fit to the distribution. The /i, a and 
\ 2 /ndf values for such distributions at centrality bins 
varying from 0-2% to 48-50% have been extracted and 
are shown in Fig. |l3[ The x 2 jndj values are seen to 
be between 1 and 2, which indicates that the distribu- 
tions are well described by Gaussians. The fluctuations 
in Et have been calculated by using Eq. 1 and are shown 
in Fig. [l4| The relative fluctuations are observed to in- 
crease in going from central to peripheral collisions. Et, 
measured by the MIRAC was used in the online trigger 
to define the most central event sample with a threshold 
which occurred in the region of the top 14% — 18% of the 
total cross section. This region is not analyzed to avoid 
trigger bias effects in the measured Et distribution. 

Et has a strong correlation with the number of par- 
ticipant nucleons or the number of effective collisions 
they undergo [[33| . In an attempt to understand the fluc- 
tuations in terms of the number of participant nucle- 
ons, the VENUS event generator has been used to de- 
termine the ratio h(Et) / fi(Npart) as a function of cen- 
trality. The /i(£t) per participant has been found to be 



~ l.lGeV ± 0.2. This is shown in Fig. ^5} The main 
sources of error here are due to the uncertainty in the 
calculation of the number of participants and the finite 
resolution of the the calorimeters. Similar results were 
also obtained from WA80 and HELIOS collaborations 
Ref. [|l]|2|]. While the WA80 Collaboration has shown 
that Et per participant is independent of projectile, tar- 
get, and centrality but depends only on the number of 
wounded nucleons and the beam energy, the WA98 Col- 
laboration has shown that transverse energy does deviate 
from a linear dependence on the number of participants 
for Pb+Pb collisions @. 

The relative fluctuations in transverse energy can be 
analyzed in a participant picture similar to that employed 
in the case of photons and charged particles [p3[]. An 
expression similar to Eq. 4 can be obtained, where the 
first term depends on the fluctuations in the transverse 
energy deposited by each particle produced per partici- 
pant nucleon, with the second term coming from impact 
parameter fluctuations within the acceptance of the de- 
tector. Since the first term depends greatly on the de- 
tector characteristics, we compare the transverse energy 
fluctuations in data to those obtained from a fast simula- 
tion of the MIRAC and ZDC characteristics in VENUS 
in which the energy resolution for each particle was ap- 
plied separately when computing the total transverse en- 
ergy [p3[]. Fig. |l](a and b) shows the comparison of the 
simulated Et and E-p distributions with those from data. 
The agreement is seen to be quite reasonable. 

Fig. [if] shows the comparison of Et fluctuations 
from data to those obtained from simulated events us- 
ing VENUS. The fluctuations are plotted as a function of 
the mean number of participants in various 2% bins of 
centrality obtained from Ep. Errors shown in the data 
are mainly due to uncertainties in the determination of 
number of participants and the finite energy resolution of 
the calorimeters as discussed earlier. It is seen that the 
fluctuation in data are systematically smaller than those 
obtained form VENUS. As discussed in Ref. [^3[| many 
effects like energy-momentum degradation of nucleonic 
objects in successive scatterings and short range correla- 
tions between nucleons in a nucleus may be responsible 
for the decrease in fluctuations in Et of data as compared 
to those obtained from simulations in VENUS. The role 
of re-scattering has also to be understood in this context. 

VII. SUMMARY 

A detailed event-by-event study of fluctuations in the 
multiplicities of charged particles and photons and trans- 
verse energy has been carried out using data from the 
WA98 experiment. This has been done varying both the 
centrality and rapidity intervals. It is observed that the 
relative fluctuations increase with increase in the impact 
parameter interval. Hence it is important to control the 
impact parameter dependence by studying narrow bins in 
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centrality. For 2% centrality bins, within which the dis- 
tributions of charged particle and photon multiplicities, 
as well as the transverse energy, are well-described by 
Gaussians, the contribution from impact parameter fluc- 
tuations is around one, as expected. The fluctuations in 
multiplicities and Et are found to increase in going from 
central towards peripheral events. A decrease in accep- 
tance has been found to result in a decreased multiplicity 
fluctuations. Using a simple statistical analysis one can 
explain the observed decrease in a smaller acceptance 
knowing the fluctuations in a larger acceptance window. 
The observed centrality dependence of the multiplicity 
fluctuations of charged particles have been found to agree 
reasonably well with results obtained from a simple par- 
ticipant model which takes into account impact param- 
eter fluctuations from VENUS and multiplicity fluctua- 
tions from NN data, within the quoted systematic er- 
rors. For photons the fluctuations are found to be slightly 
lower compared to those obtained from the participant 
model. 

Similar calculations have been performed for the 
transverse energy and multiplicity distributions using 
VENUS. The transverse energy fluctuations from exper- 
imental data are found to be smaller than those observed 
in VENUS. On the other hand, after corrections for 
charged particle contamination in the photon-like clus- 
ters, the relative fluctuations of photons appear to be in 
rather good agreement with VENUS. 

Acknowledgements 

We wish to express our gratitude to the CERN accel- 
erator division for the excellent performance of the SPS 
accelerator complex. We acknowledge with appreciation 
the effort of all engineers, technicians and support staff 
who have participated in the construction of this experi- 
ment. 

This work was supported jointly by the German 
BMBF and DFG, the U.S. DOE, the Swedish NFR and 
FRN, the Dutch Stichting FOM, the Polish KBN under 
Contract No. 621/E-78/SPUB-M/CERN/P-03/DZ21 1/, 
the Grant Agency of the Czech Republic under contract 
No. 202/95/0217, the Department of Atomic Energy, 
the Department of Science and Technology, the Coun- 
cil of Scientific and Industrial Research and the Uni- 
versity Grants Commission of the Government of India, 
the Indo-FRG Exchange Program, the PPE division of 
CERN, the Swiss National Fund, the INTAS under Con- 
tract INTAS-97-0158, ORISE, Grant-in-Aid for Scien- 
tific Research (Specially Promoted Research & Interna- 
tional Scientific Research) of the Ministry of Education, 
Science and Culture, the University of Tsukuba Special 
Research Projects, and the JSPS Research Fellowships 
for Young Scientists. ORNL is managed by UT-Battelle, 
LLC, for the U.S. Department of Energy under contract 
DE-AC05-00OR22725. The MIT group has been sup- 
ported by the US Dept. of Energy under the cooperative 
agreement DE-FC02-94ER40818. 



[1] L. van Hove, Phys. Lett. B118, 138 (1982); H. Heisel- 
berg, Phys. Rep. 351, 161 (2001); E.V. Shuryak, Phys. 
Lett. B423, 9 (1998); M. Gazdzicki and S. Mrowczynski, 
Z. Phys. CS7, 127 (1992); S. Mrowczynski, Phys. Rev. C 
57, 1518 (1998); Phys. Lett. B430 9 (1998). 

[2] Proceedings of Quark Matter '99, Nucl. Phys A 661 
(1999). 

[3] Proceedings of Quark Matter '2001, To appear in Nucl. 
Phys A. 

[4] M. Stephanov, K. Rajagopal and E. Shuryak, Phys. Rev. 

Lett. 81, 4816 (1998); Phys. Rev. D 61, 114028 (1999). 
[5] G. Baym and H. Heiselberg, Phys. Lett. B469, 7 (1999). 
[6] MAsakawa, U.Heinz, and B.Miiller, Phys. Rev. Lett. 85, 

2072 (2000). 

[7] S. Jeon and V. Koch Phys. Rev. Lett. 83, 5435 (1999); S. 

Jeon and V. Koch, Phys. Rev. Lett. 85, 2076 (2000). 
[8] L. van Hove, Z. Phys. C21, 93 (1984); J.I. Kapusta and 

A.P. Vischer, Phys. Rev. C52, 2725 (1995). 
[9] K. Rajagopal and F. Wilczek, Nucl. Phys. B399, 395 
(1993); Nucl. Phys. B404, 577 (1993); J.D. Bjorken, Int. 
J. Mod. Phys. A7, 4189 (1992); T.C. Brooks et al., (Mini- 
max Collaboration), Phys. Rev. D61, 032003 (2000); J. -P. 
Blaizot and A. Krzywicki, Phys. Rev. D46, 246 (1992). 
[10] M.M. Aggarwal et al., (WA98 Collaboration) Phys. Rev. 
C64, 011901R (2001). 
T. Akesson et al., Z. Phys. C38, 383 (1988). 
T. Akesson et al., Phys. Lett. B214, 214 (1988). 
T. Alber et al., (NA49 Collaboration) Phys. Rev. Lett. 75, 
3814 (1995). 

[14] H. Appelshauser et al., (NA49 Collaboration), Phys. Lett. 
B459, 679 (1999); G. Roland, Nucl. Phys. A638, 91c 
(1998). 

[15] G. Baym et al., Phys. Lett. B219, 205 (1989). 

[16] F. Corriveau et al., (NA34 Collaboration), Z. Phys. C38, 
15 (1988); J. Schukraft et al., (NA34 Collaboration), Z. 
Phys. C38, 59 (1988). 

[17] R. Albrecht et al., Z. Phys. C45, 31 (1989). 

[18] Proposal for a Large Acceptance Hadron and Pho- 
ton Spectrometer, Preprint CERN/SPSLC 91-17, SP- 
SLC/P260. 

[19] M.M. Aggarwal et al., (WA98 Collaboration), Phys. Lett. 

B420, 169 (1998). 
[20] W.T. Lin et al., Nucl. Instrum. Meth. A389, 415 (1997). 
[21] M.M. Aggarwal et al., Nucl. Instrum. and Methods in 

Phys. Res. A424, 395 (1999). 
[22] M.M. Aggarwal et al., (WA98 Collaboration), Phys. Lett. 

B458, 422 (1999). 
[23] T. C. Awes et al., Nucl. Instrum. Methods Phys. Res. Sect. 

A 279, 497 (1989). 
[24] K. Werner, Phys. Rep. C232, 87 (1993). 
[25] D.P. Mahapatra, B. Mohanty and S.C. Phatak, 



[11] 
[12] 
[13] 



print:nucl-ex/010801 1 ; To be published in Int. J. Mod. 



Phys. A. 

[26] Rudolph C. Hwa, Phys. Lett. B201, 165 (1988). 



9 



[27] A. Giovannini and L. Van. Hove, Z. Phys. C30, 391 
(1986). 

[28] M.M. Aggarwal et al., (WA98 Collaboration), Eur. Phys 

J. C18, 651 (2001). 
[29] G.J. Alner et al., (UA5 Collaboration), Phys. Rep. C154, 

247 (1987). 

[30] J. Whitmore, Phys. Rep. C27, 188 (1976); J. Whitmore, 

Phys. Rep. CIO, 274 (1974). 
[31] B. Nilsson-Almqvist, E. Stenlund, Comp. Phys. Comm. 

43 (1987) 387. 

[32] M.M. Aggarwal et al., Nucl. Instrum. and Methods in 

Phys. Res. A372, 143 (1996). 
[33] Gordon Baym et al., Phys. Lett. B190, 29 (1987); A.D. 

Jackson and H. Bogglid, Nucl. Phys. A470, 669 (1987). 



10 





a riiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii.il 
50 100 150 200 250 300 350 400 450 500 

E T (GeV) 

FIG. 1. Minimum bias distributions of (a) transverse energy Et, (b) forward energy Ef produced in Pb induced reactions at 
158-A GeV on Pb. Solid histograms show the results obtained from VENUS event generator, (c) shows the anti-correlation of 
measured total transverse energy, Et, and forward energy, Ep. 
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FIG. 2. Minimum bias distributions of (a) 7-like cluster multiplicity, and (b) charged particle multiplicity produced in Pb 
induced reactions at 158-^4 GeV on Pb. The multiplicity distributions for the top 1%, 2%, and 5% most central events are also 
shown and fitted to Gaussian distributions. 
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FIG. 3. Gaussian fit parameters of the multiplicity distributions of 7-like clusters and charged particles for increasing centrality 
bins of increasing width. The centrality selection has been made by increasing the widths of Et bins corresponding to — 1%, 
— 2%, — 3%,- • -,0 — 10% of the minimum bias cross section. 



600 
^ 500 
^4 400 
^ 300 
^ 200 
zL 100 




u 



IN 







»*..« 



• • « • 



0-2 8-10 18-20 28-30 38-40 48-50 58-60 

% cs ( in 2 % bins ) 





600 




500 




400 


o 


300 


ZL 200 




too 









40 




30 


o 






20 




10 











o 






3 


T3 






2 


M 









••?.. 

I I I I I 



♦♦♦♦♦♦ 



♦♦♦♦ 



0-2 8-10 18-20 28-30 38-40 48-50 58-60 

% cs ( in 2 % bins ) 



FIG. 4. Gaussian fit parameters of the multiplicity distributions of 7 — like clusters and charged particles as a function of 
centrality. The centrality selection has been made by selecting 2% bins in minimum bias cross section, viz., — 2%, 2 — 4%, 
4 — 6%,- • -,62 — 64%. The multiplicity distributions for these centrality bins are near perfect Gaussians as can be seen from the 
X 2 /ndf values. 
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FIG. 5. Fluctuations of the multiplicity of 7-like clusters and charged particles within the full coverage of PMD and SPMD for 
various 2% bins of the minimum bias cross section. 
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FIG. 6. Centrality dependence of the Gaussian fit parameters of the multiplicity distribution of 7-like clusters and charged 
particles within the common coverage of PMD and SPMD. 
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FIG. 7. Centrality dependence of the fluctuations of the multiplicity of photons and charged particles within the common 
coverage of the PMD and SPMD. 
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FIG. 8. Multiplicity fluctuations of photons and charged particles for two r\ acceptance selections. The open squares represent 
estimated values fluctuations in 0.5 unit of 8q from the observed fluctuations in 1.0 unit of 5r). 
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FIG. 9. Variation of //, a, and \ 2 /ndf of the distribution of the number of participants as a function of centrality. 
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FIG. 10. The relative fluctuations, uj c h, of the charged particle multiplicity as a function of number of participants, 
experimental data are compared to calculations from a participant model and those from VENUS event generator. 
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FIG. 11. The relative fluctuations, lu 1 of photons as a function of number of participants. The data presented show the 
fluctuations in 7 — like clusters and photons after correction for efficiency and purity. These are compared to calculations from a 
participant model and those from VENUS event generator. 
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FIG. 13. Centrality dependence of fx, a, and \ 2 /ndf of the transverse energy distribution. The centrality selection is base on 
E-p measured with the ZDC. 
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FIG. 15. Et per participant as a function of centrality. The vertical solid line indicates the estimated systematic error in E*r 
per participant. 
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FIG. 16. Centrality dependence of the relative fluctuations in transverse energy, Et, with the centrality selected from the Ef. 
The result is compared to VENUS using similar centrality selection criteria. 
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